Novel NN interaction and the spectroscopy of light nuclei 
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Nucleon-nucleon (NN) phase shifts and the spectroscopy of A < 6 nuclei are successfully described 
by an inverse scattering potential that is separable with oscillator form factors. 
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Nucleon-nucleon (NN) potentials that describe avail- 
able two-body data have a long and multi-faceted his- 
tory. High precision fits have improved with time even 
as more precise experimental data have become avail- 
able. Three-nucleon (NNN) potentials have a shorter 
history but are intensively investigated at the present 
time. Disparate foundations for these potentials, both 
NN & NNN, have emerged. On the one hand, one sees 
the predominant meson-exchange potentials sometimes 
supplemented with phenomcnological terms to achieve 
high accuracy in fitting NN data (Bonn 0, Nijmegen 
0, Argonne 0, Idaho 0, IS 0) and NNN data (Ur- 
bana 0,0, Illinois 0], Tucson-Melbourne 0E3)- 0n 
the other hand, one sees the emergence of potentials with 
ties to QCD which are either meson- free ^lj; or inter- 
twined with meson-exchange theory 0, . 

All these potentials are being used, with unprecedented 
success, to explain a vast amount of data on light nuclei 
in Quantum Monte Carlo ap proa ches and ab initio 
no-core shell model (NCSM) [H The overwhelm- 
ing success of these efforts have led some to character- 
ize these approaches as leading to a 'Standard Model' of 
non-relativistic nuclear physics. 

Chief among the outstanding challenges is the compu- 
tational intensity of using these NN + NNN potentials 
within the presently available many-body methods. For 
this reason, most ab initio investigations have been lim- 
ited to A < 12. The situation would be dramatically 
simpler if cither the NN potential alone would be suffi- 
cient or the potentials would couple less strongly between 
the low momentum and the high momentum degrees of 
freedom. If both simplifications are obtained, the future 
for applications is far more promising. 

In the present work, we derive and apply a new class of 
potentials which have no apparent connection with the 
two well-established lines of endeavor. We develop J- 
matrix inverse scattering potentials (JISP) that describe 
NN data to high accuracy and, with the off-shell free- 
dom that remains, we obtain excellent fits to the bound 
and resonance states of light nuclei up to A = 6. Our 
NN off-shell freedom is sufficient to describe these lim- 
ited data without the need for NNN potentials. As an 
important side benefit, we find that these potentials lead 
to rapid convergence in the ab initio NCSM evaluations 
presented here. We hope that these potentials will open 



a fruitful path for evaluating heavier systems and spur 
the development of extensions to scattering problems. 

It is important to stress that our NN potentials have 
the same symmetries as the conventional NN potentials 
mentioned above (without charge symmetry breaking at 
present), but they are not constrained by meson exchange 
theory, by QCD or by locality. This does not mean our 
NN potentials are inconsistent with those constraints, 
however. 

By means of the J-matrix inverse scattering approach 
(lif we construct NN potentials as matrices in an oscil- 
lator basis with Hui — 40 MeV using the Nijmegen np 
phase shifts [Tt^ . Following Ref. [HI, we obtain inverse 
scattering tridiagonal potentials (ISTP) that are tridi- 
agonal (quasi-tridiagonal) in uncoupled (coupled) partial 
waves. The dimension of the potential matrix is specified 
by the maximum value of N — In + I and is referred to 
as an Nfku potential. In order to improve the descrip- 
tion of the phase shifts, we develop a 9fio;-ISTP in odd 
waves instead of the 7fio;-ISTP of Ref. ^^|- We retain 
an 8fitJ-ISTP in the even partial waves. 

Next we perform various phase equivalent transforma- 
tions (PETs) of the obtained ISTP. In the coupled sd 
waves, we perform the same PET as in Ref. |15| but 
with different rotation angle i? = 11.3° to improve the 
description of the deuteron quadrupole moment Q. We 
then find improvement in 3 H and He binding energies. 
We also perform similar PETs mixing lowest oscillator 
basis states in the 3 p2, 3 pi, 3 d-2 and 1 p\ waves with the 
rotation angles of § — +8°, —6°, +25° and —16° respec- 
tively to improve the description of the 6 Li spectrum. 
The obtained interaction fitted to the spectrum of A = 6 
nuclei, is refered to as JISP6. The non-zero matrix ele- 
ments of the JISP6 interaction are presented in TablesQ]- 
13 (in hw = 40 MeV units). 

The deuteron properties provided by JISP6 are com- 
pared with those of some other realistic potentials in Ta- 
ble EH 

We perform calculations of light nuclei in the NCSM 
with JISP6 plus the Coulomb interaction between pro- 
tons. To improve the convergence, we perform the Lee- 
Suzuki transformation to obtain a two-body effective in- 
teraction as is discussed in Ref. jj^j. We obtain the 
effective interaction in a new basis (hoj — 15 MeV) 
within an N max Hu) model space where N max signifies 
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TABLE I: Non-zero matrix elements in hu = 40 MeV units TABLE V: Non-zero matrix elements in huj — 40 MeV units 
of the JISP6 matrix in the x sq partial wave. of the JISP6 matrix in the po partial wave. 







y( _ yl 

'n,n + l *n + l,n 


n 


V 1 

y nn 


yi _ yl 

'n,n + l *n + l,n 


-0.370829835410 




0.132663053234 





-0.143164548564 


0.020755069079 


1 -0.148826473896 




0.006448104375 


1 


0.082988173615 


-0.120094506156 


2 0.152835073179 




-0.120193538281 


2 


0.310447079465 


-0.116102071897 


3 0.187138532148 




-0.029504403821 


3 


0.065044984944 


0.013609203897 


4 -0.005584124194 






4 


-0.026555044004 




TABLE II: Non-zero matrix elements 


in 


huj = 40 MeV units 


TABLE VI: Non-zero matrix elements in 


huj = 40 MeV units 


of the JISP6 matrix in the 1 pi partial 


wave. 


of the JISP6 matrix m the pi partial wave. 


n Vnn Vn, n + 1 — V n +1, n 


Vn,n + 2 — Kl+2, n 


n 


V 1 V 1 .-, = V' ; ,i 

v nn ' 7i, n + 1 * rt-\- 1 , n 


Vn,n+2 — Ki+2,n 


0.631081576474 -0.251382936855 


0.413319237868 





0.249679784904 -0.164761352606 


0.157602869176 


1 -0.293390247307 -0.118539824524 




1 


0.044327922667 -0.176615480839 




2 0.454133632867 -0.230186013455 




2 


0.514099248363 -0.275733929941 




3 0.348035837559 -0.090043227032 




3 


0.423324941377 -0.108223480414 




4 0.049221178231 






4 


0.055397268090 




TABLE III: Non-zero matrix elements 


in 


ftw = 40 MeV units 


TABLE VII: Non-zero matrix elements in 


ftw = 40 MeV units 


of the JISP6 matrix in the 1 d,2 partial 


wave. 


of the JISP6 matrix m the d,2 partial wave. 


n V l nn 




V 1 — V 1 

" n, n-\-± v n+l,n 


n 


V 1 V 1 — V 1 

v nn v n,n-l-l v n + l,n 




-0.040699390045 




0.037531685348 





-0.662113235725 0.759732268972 


-0.571851583929 


1 -0.111761745807 




0.069791608541 


1 


0.175448232529 0.273660320760 




2 -0.134996618188 




0.045265020638 


2 


-0.263880156721 0.086029122705 




3 -0.031270631267 






3 


-0.063231092747 




TABLE IV: Non-zero matrix elements 


in 


hu) = 40 MeV units 


TABLE VIII: Non-zero matrix elements 


in hw = 40 MeV 


of the JISP6 matrix in the 1 f$ partial 


wave. 


units of the JISP6 matrix in the 3 /3 partial wave. 






V 1 — V 1 

"n, n + 1 * 71+1, 71 


n 


V 1 

y nn 


V 1 — V 1 

'n, n + 1 *n + l,n 


0.019468932284 




-0.018631285428 





0.026326206898 


-0.014285757490 


1 0.083569595544 




-0.055488297878 


1 


0.035674429367 


-0.016797566427 


2 0.106821875609 




-0.032273326919 


2 


0.028543592124 


-0.008290586001 


3 0.021063860199 






3 


0.006036946613 





the many-body oscillator basis cutoff. The results of our 
NCSM calculations for binding energies of 3 H, 3 He (in the 
lifiw model space), 4 He (in the 12hui model space), 6 He 
(in the 8huj model space) and 6 Li (in the lOftw model 
space) nuclei are compared in Table IXTTl with the calcu- 
lations in various approaches [Faddeev, Green's-function 
Monte Carlo (GFMC), NCSM] with realistic NN [CD- 
Bonn, Nijmegen-I (Nijml), Nijmegen-II (Nijmll), and Ar- 
gonne (AV18 and AV8')] and NNN [Urbana (UrblX) 
and Tucson-Melbourne (TM and TM')] potentials. To 
give an estimate of the convergence of our calculations, 
we present the difference between the given result and 
the result obtained in the next smaller model space in 
parenthesis after our JISP6 results. It is seen that the 
convergence of our calculations is adequate. 



The convergence patterns are also illustrated by Fig.Q] 
where we present the huj dependence of the 6 Li ground 
state energy in comparison with the results of Ref. [lflj 
obtained in NCSM with CD-Bonn interaction. The hw 
dependence with the JISP6 interaction is weaker over a 
wide interval of hu> values. This is a signal that conver- 
gence is improved relative to CD-Bonn. The variational 
principle cannot be applied to the NCSM calculations 
with effective interactions so the convergence may be ei- 
ther from above or below. However, we may surmise that 
the residual contributions of neglected three-body effec- 
tive interactions is more significant in the CD-Bonn case. 

Returning to the results presented in Tables IXIHXIII 
we see that the JISP6 interaction provides a realistic de- 
scription of the ground states of light nuclei competitive 
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TABLE IX: Non-zero matrix elements in hw = 40 MeV units 
of the JISP6 matrix in the sd coupled waves. 



V^i matrix elements 



v n 



v, 



— y-n+i 



-0.508274040822 

1 -0.276168029473 

2 -0.009473803659 

3 0.152873734289 

4 0.037547929880 

ydc 

n V* 



0.214156446587 
0.080907735691 
-0.051881443108 
-0.055193589842 



dd 

n,n-\-l 



-trdd 

' ' *n-\-l,n 



0.050878349132 
0.322126471805 
0.308516673061 
0.061200037193 



-0.094173649477 
-0.178808793641 
-0.093012604766 



1 n. 



V, 
V„ d 



sd 



v, 



sd 



V, 



ds 



v n,n + l 



-0.048516342749 
0.068044496963 
0.049400578816 

-0.001503998139 



-0.411771355351 
-0.060476585693 
-0.080187106458 
-0.020205646231 



0.205731982741 



TABLE X: Non-zero matrix elements in hw = 40 MeV units 
of the JISP6 matrix in the pf coupled partial waves. 



V vv , matrix elements 



V 



v- 



n,n+l 



7i-\-L,n 



ypp _ ypp 



-0.257052769018 
0.035950531483 
-0.209221226255 
-0.151546344031 
-0.015878201100 

rff 



0.215269922150 
0.103784437050 
0.103221679634 
0.037329967115 



-0.171332097424 



V i matrix elements 



yff 



yff _ yff 



-0.019836117401 

1 -0.010058323807 

2 0.001646202464 

3 0.000388504318 



0.008292672214 
0.000628665253 
-0.000973797731 



ypf 



V 



fp 



V , matrix elements 

nn 

ypf — yfp 

n v nn v n.n. 



v pf ri = v f L 





1 -0.026186898553 

2 -0.024757588981 

3 -0.014708906826 

4 0.000024653107 



0.018163699530 
0.023478346345 
0.023707438623 
0.006271279847 



0.003299666393 
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FIG. 1: (Color online) hto dependence of the 6 Li ground state 
energy obtained with JISP6 interaction in comparison with 
the one obtained in NCSM with CD-Bonn potential |Tsft . 



with the quality of descriptions previously achieved with 
both NN and NNN forces. 

This conclusion is supported by the calculations of the 
spectra of A = 6 nuclei with Hui = 15 MeV presented 
in Table IXlTTl We again present in parenthesis the dif- 
ference between the given excitation energy and the re- 
sult obtained in the next smaller model space. It is seen 
that the 6 Li spectrum is well-reproduced in our calcu- 
lations. The most important difference with the experi- 
ment is the excitation energy of the (1^,0) state. How- 
ever £^(1^,0) goes down rapidly when the model space 
is increased and better results are anticipated in a larger 
model space. The JISP6 results for 6 Li spectrum are also 
seen to be competitive with results from modern realistic 
NN + NNN interaction models. We note here that the 
Li spectrum was found to be significantly sensitive 
to the presence of the NNN force and this motivated our 
adoption of 6 Li for these comparisons. 

We return to the underlying rationale for our approach 
and ask why it is conceivable that an NN interaction 
alone may be as successful as the NN + NNN potentials 
mentioned at the outset. That this is feasible may be ap- 
preciated from the theorem of Polyzou and Glockle [20( . 
They have shown that changing the off-shell properties 
of two-body potentials is equivalent to adding many- 
body interactions. This theorem coupled with our lim- 
ited results suggests that our inverse scattering NN po- 
tential plus off-shell modifications is roughly equivalent, 
for the observables so far investigated, to the successful 
NN + NNN potential models. 

Clearly, more work will be needed to carry this to nu- 
clei with A>7 and see if the trend continues. Based on 
the results presented, the additional off-shell freedoms 
remaining may well serve to continue this line of fit- 
ting properties for some time. When it eventually breaks 
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TABLE XI: JISP6 deuteron property predictions in comparison with the ones obtained with various realistic potentials. 



Potential 



E d , MeV 



d state 
probability, % 



rms radius, 
fm 



Q, fm 2 



As. norm, 
const, s/ s , 
fm" 1 / 2 



V 



s>f» 



JISP6 
Nijmegen-II 

AV18 
CD-Bonn 
Nature 



-2.224575 
-2.224575 
-2.224575 
-2.224575 
-2.224575(9) 



4.1360 
5.635 
5.76 
4.85 



1.9647 
1.968 
1.967 
1.966 
1.971(6) 



0.2915 
0.2707 
0.270 
0.270 
0.2859(3) 



0.8629 
0.8845 
0.8850 
0.8846 
0.8846(9) 



0.0252 
0.0252 
0.0250 
0.0256 
0.0256(4) 



TABLE XII: The binding energies of 3 H, 3 He, 4 He, 6 He and 6 Li nuclei. 



Potential 


3 H 


3 He 


4 He 


6 He 


6 Li 


JISP6, NCSM 


8.461(5) 


7.751(3) 


28.611(41) 


29.072(69) 


31.48(27) 


CD-Bonn+TM, Faddeev [17] 


8.480 


7.734 


29.15 






AV18+TM, Faddeev [17] 


8.476 


7.756 


28.84 






AV18+TM, Faddeev [17] 


8.444 


7.728 


28.36 






Nijml+TM, Faddeev [171 


8.392 


7.720 


28.60 






Nijmll+TM, Faddeev [17] 


8.386 


7.720 


28.54 






AV18+UrbIX, Faddeev [17] 


8.478 


7.760 


28.50 






AV18+UrbIX, GFMC [7] 


8.47(1) 




28.30(2) 


27.64(14) 


31.25(11) 


AV8'+TM', NCSM [18] 








28.189 


31.036 


Nature 


8.48 


7.72 


28.30 


29.269 


31.995 



TABLE XIII: Excitation energies E x (in MeV) of A = 6 nu- 
clei. 



6 Li 


Nature 


JISP6 


AV8'+TM' 


AV18+UrbIX 


Model space 




whuj 


NCSM, 6hw [18] 


GFMC [7] 




0.0 


0.0 


0.0 


0.0 


E x {3+,0) 


2.186 


2.102(4) 


2.471 


2.72(36) 


^(0+1) 


3.563 


3.348(24) 


3.886 


3.94(23) 


E x (2+,0) 


4.312 


4.642(2) 


5.010 


4.43(39) 


E x (2+,1) 


5.366 


5.820(4) 


6.482 




E.(lt,0) 


5.65 


6.86(36) 


7.621 




6 He 


Nature 


JISP6 


AV8'+TM' 


AV18+UrbIX 


Model space 




8huj 


NCSM, Qhto [18J 


GFMC [JJ 


E x (0+,1) 


0.0 


0.0 


0.0 


0.0 


E x (2+,1) 


1.8 


2.59(13) 


2.598 


1.80(18) 



down, NNN potentials may be needed. 
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